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Abstract 

We employ the language of Cartan's geometry to present a model 
for studying vector spaces of Killing two-tensors denned in pseudo- 
Riemannian spaces of constant curvature under the action of the cor- 
responding isometry group. We also discuss geometric properties of 
joint invariants of Killing two-tensors defined in the Euclidean plane 
to formulate and prove an analogue of the Weyl theorem on joint in- 
variants. In addition, it is shown how the joint invariants manifest 
themselves in the theory of superintegrable Hamiltonian systems. 



1 Introduction 

The purpose of this article is three- fold. First, we wish to introduce a general 
model which forms a suitable framework for describing the invariant theory 
of Killing tensors from the viewpoint of Cartan's geometry (see [T2l [131 [T5| 
35] for more details), in particular its moving frame method and the 
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theory of fiber bundles. Apparently, the idea to study the vector spaces of 
second order differential operators defined in the Euclidean plane under the 
action of the corresponding isometry group goes back to the classical 1965 
Lie group theory paper by Winternitz and Fris [13]. In 2002 an analogous 
idea resurfaced in the problem of studying the isometry group action in a 
vector space of Killing two-tensors in the paper by McLenaghan, Smirnov 
and The [2B] on classical Hamiltonian systems. A series of recent works on 
the subject represent a steady development of the invariant theory of Killing 

tensors PQ El dSl HZl I2SI I2ZI I2H1 IS01 ISU EZl I3SI ISl SS]- 

Our second goal is to extend the study of joint invariants of Killing ten- 
sors introduced in Smirnov and Yue [35] • More specifically, we formulate and 
prove an analogue of the Weyl theorem on joint invariants in classical invari- 
ant theory (see Olver [33] and Weyl [12] for more details). Furthermore, we 
introduce the concept of a resultant of Killing tensors, which is also inspired 
by its analogue in classical invariant theory (see Olver [33] for more details). 

Finally, we employ the resultants to derive a new characterization of 
superintegrable systems. Recall that although examples of superintegrable 
systems have been known since the time of Kepler, a systematic development 
of a general theory of superintegrable systems originated in the pioneering 
works by Winternitz and collaborators JHJ [25] (see also Smirnov and Winter- 
nitz [39J for more references). The theory has been actively developed since 
then for both quantum and classical Hamiltonian systems. Thus, a general 
classification and structure theory for superintegrable systems defined in two- 
and three-dimensional conformally flat spaces has been introduced in a series 
of recent papers by Kalnins, Kress and Miller (see [19], [20l [21], [221 [23] and 
the relevant references therein). As is well-known, a wide range of methods 
are employed to study superintegrable systems. We shall demonstrate that 
the geometric methods stemming from the invariant theory of Killing tensors 
also prove their worth in this study and present a new prospective. Thus, 
the main theorem presented in this work is a geometric characterization of 
the Kepler potential in terms of the existence of a vanishing resultant of the 
associated Killing two-tensors. 

2 The model 

In what follows, we shall assume that (M, g) is an m-dimensional (pseudo- 
) Riemannian manifold of constant curvature. Generalized Killing tensors 
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can be denned as the elements of the vector space of solutions to the fol- 
lowing overdetermined system of PDEs, called the generalized Killing tensor 
equation, given by 

[[...[K,g],g],...,g] = (n+1 brackets), (1) 

where [ , ] denotes the Schouten bracket [31] and K is a symmetric con- 
travariant tensor of valence p. The generalized Killing tensor equation ([1]) is 
determined by given m > 1, n > and p > 0. Let n, m and p be fixed and 
K^{M) denote the corresponding vector space of solutions to the generalized 
Killing tensor equation (JTJ). The Nikitin-Prylypko-Eastwood (NPE) formula 

PE1E2] 

d = dunJC?(M) = ^±±( P + m - 1 )( P + n + m ), (2) 
m \ m — 1 / \ m — 1 / 

represents the dimension d of the vector space /C^(M). The NPE formula (J2J) 
generalizes the Delong-Takeuchi-Thompson (DTT) formula [31101111] derived 
for the case n — 0. When (M, g) is not of constant curvature, d given by 
d2]) is the least upper bound for the dimension of the space of solutions to 
(PQ). Of particular importance, due to their relevance in various problems of 
mathematical physics, are the elements of the vector space /Cq(M) ~ M. d for 
a given space (M, g) (see, for example, [H [3 UJil [HI [221 EHl EHl S3] 

and the relevant references therein). 

The geometric study of the elements of /Cq(M) having distinct (and real 
in the case of (M, g) being pseudo-Riemannian) eigenvalues and hypersurface 
forming eigenvectors (eigenforms) is pivotal in the Hamilton- Jacobi theory 
of orthogonal separation of variables (see [21 [El [13, [221 [231 [2S1 EE] and the 
relevant references therein). An element K e JCq(M) enjoying these proper- 
ties generates an orthogonal coordinate web (see JTOj [21 E] and the relevant 
references therein) which consists of m foliations, the leaves of which are 
m — 1-dimensional hypersurfaces orthogonal to the eigenvectors of K. In 
this way the orthogonal coordinate web is adapted to the hypersurface form- 
ing eigenvectors (eigenforms) of the corresponding Killing tensor K. The 
most natural framework for such a study is best described in terms of an ap- 
propriate Cartan's (or Klein's) geometry [351 (HI EH] > depending on a given 
(pseudo-) Riemannian manifold (M, g) and its geometry. Indeed, let G be 
the isometry group of (M, g), and H C G - a closed subgroup of the Lie group 
G. The study of the orthogonal coordinate webs generated by the elements 
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of /Cq(M) with the properties prescribed above, is established in the study 
of the principal fiber if -bundle: 7i"i : G — > G/H ~ M, where the left coset 
space G/H is identified with M on which G acts transitively. We also have 
the natural structure of a vector bundle: 7r 2 : Kq(M) — > G/H ~ M. The 
orthogonal coordinate webs are thus defined in the homogeneous space M 
and the basic equivalence problem can be formulated as follows: Given two 
orthogonal coordinate webs Wi, W2 G M generated by the corresponding el- 
ements Ki,K 2 G /Cg(M) respectively. The coordinate webs Wi, W2 are said 
to be equivalent iff there exists a group element g & G and £ G K such that 
pKi = K 2 +£g. Thus, given two orthogonal coordinate webs Wi, i = 1,2, one 
wishes to know whether or not they are equivalent as defined above. More 
specifically, the transitive action G O M induces the corresponding non- 
transitive action G O Kq(M), solving the equivalence problem then amounts 
to the study of the orbit space Kq(M)/G. This structure leads to the princi- 
pal G-bundle tt 3 : £§(M) -> K%(M)/G, provided the action G O Kg(M) is 
free (for more details on the fiber bundle theory see, for example, Fatibene 
and Francaviglia [H]). Finally, one can define a map / : /Cq(M)/G — > G, so 
that the following diagram commutes. 



G. 



11 ^G/H ~ M 



./ 



£g(M)/G 



713 



7T 2 



Now let x E G/H and an element K G /C§(M) has pointwise distinct eigen- 
values and hypersurface forming eigenvectors. If x is not a singular point 
of the coordinate web (a point where the eigenvalues of K coincide), we can 
use the metric g to orthonormalize the eigenvectors Ex, ... , E m of K and 
then use this frame of eigenvectors as a basis of T^G/H, which has a group 
conjugate to H acting on it. Thus, each fiber tt ] ~ 1 (x) can be identified with 
an orthonormal frame of eigenvectors of K. In the other direction the fiber 
7r 2 _1 (x) can be identified with the corresponding vector space in /Cq(M) for 
x G M. Furthermore, in this vector space we can fix the parameters so that 
the resulting element of the vector space is precisely the Killing tensor K 
above evaluated at x G M. It is now a point in the fiber 7r 2 _1 (x). Next, 
under the projection 7^3 this point is mapped to the corresponding point in 
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the orbit space /Cq(M)/G, which is nothing, but the orbit generated by K 
under the action of G O /Cq(M). Finally, choosing a function / lifting the 
action to G is equivalent to either choosing a cross-section K through the 
orbits of K.q(M)/G or fixing a frame. Then (local) invariants of the group 
action G O JCq(M) are the coordinates of the canonical forms obtained as 
the intersection of the orbits with an appropriately chosen cross-section K 
[53] . In the latter case one can solve the equivalence problem by employing 
Cartan's language of the exterior differential forms [T5J [35J, HJ [18] , while in the 
former case the problem can be solved by using the techniques of the geomet- 
ric version of the classical moving frames method as developed by Fels and 
Olver [TSJ EH [33 J (for more details of how these two approaches interact see 
|38j). The equivalence of the two approaches is manifested by the fact that 
the group G acts transitively on the bundle of frames. We note however, that 
in practice the most effective approach is to use both methods. Thus, the 
former method can be used to find canonical forms of the orbits of JCq(M)/G, 
while the latter approach comes into play when, for example, one needs to 
determine the moving frames map for a given representative of an orbit of 
K.q(M)/G. Note also that the composition map 7 = tt 3 o / : JCq(M) — ► G 
is the moving frame [T5l [T2| [T31 [33] , corresponding to the cross-section (or 
the frame of eigenvectors) prescribed by a chosen map / : KLq{M)/G — > G. 
Since "The art of doing mathematics consists in finding that special case 
which contains all the germs of generality." (D. Hilbert), in the example 
that follows we shall do just that. 

3 Illustrative example 

Consider the following example. M = E 2 , K$(M) = /C 2 (E 2 ), G = SE(2) 
and H = SO(2). This example has been extensively studied in the litera- 
ture from different points of view [U [231 [211 [36], [38] 03] . For the first time 
the orbit problem SE(2) O /Cq(E 2 ) was considered by Winternitz and Fris 
[13] in the context of a classification of quadratic differential operators which 
commute with the Laplace operator. We consider the problem of classifica- 
tion of the orthogonal coordinate webs generated by non-trivial elements of 
the vector space /Cq(E 2 ) in the framework of the above diagram. Note first 
that in this Killing tensor K G /Cq(E 2 ) with pointwise distinct eigen- 

values has necessarily hyper surf ace-forming eigenvectors, which generate one 
of the four orthogonal coordinate webs, namely elliptic-hyperbolic, polar, 
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parabolic and cartesian. The latter depends on the corresponding orbit of 
K%(E 2 )/SE(2) that K belongs to. More specifically, we fix at x G E 2 the 
following orthonormal frame fa; E\, E 2 ), where Ei,E 2 are the eigenvectors 
of K. We also assume that x G E 2 is not a singular point of K. In this 
view, x = 7Ti : SE(2) — > SE{2) / SO{2). The vectors E\,E 2 form a basis of 
T X E 2 . Let E 1 , E 2 be the corresponding dual basis of T^E 2 . Then E x , E 2 are 
horizontal in the fibration and in this frame the components of the metric g 
of E 2 and K are given by 

g ab = 5 ab E a QE b and K ab = \ a 5 ab E a & E b , a, b =1,2 (3) 

respectively, where is the symmetric tensor product, 5 ab is the Kronecker 
delta and Ai, A2 are the eigenvalues of K. Note that we have used the metric 
to orthonormalize the frame {E\,E 2 ), while the formulas ([3]) represent the 
pull-back under x = m of the metric g and the Killing tensor K to the 
bundle of frames. Moreover, we also note that by fixing a frame (Ei,E 2 ) 
we have fixed the corresponding map / : /Cq(E 2 )/ SE(2) — > SE(2). Then 
the problem of classification of the orthogonal coordinate webs generated by 
the non-trivial elements /C 2 (E 2 ) can be solved by making use of the Cartan 
theory of exterior differential forms [H [38] . Thus, upon introduction of the 
connection coefficients T of the Levi-Civita connection V via: 

V Ea E b = T ab c E c , V E E b = TjE d , a, b,c,d= 1, 2, (4) 

we arrive at the Cartan structure equations 

dE a + oA A E b = T a = 0, du a b + u a c A u c b = Vl a b = 0, (5) 

where u a b := T cb a E c are the connection one- forms, T a = \{r a bc E b A E a is the 
vanishing (V is Levi-Civita) torsion two-form and fl a b := \R a bc dE c A E d is 
the vanishing (in view of flatness of E 2 ) curvature two- form. The choice of 
the connection is not arbitrary. As is well-known from Riemannian geometry, 
given a connection V on a manifold M one can parallel propagate frames. 
For any path r between two points of M, parallel transport along r defines 
a linear mapping L(t) between the tangent spaces of two points. This linear 
map is an isometry if the connection V is a Levi-Civita connection. Clearly, 
the linear map L(r) induced by a Levi-Civita connection V maps orthonor- 
mal frames to orthonormal frames. The conditions E a A dE a = 0, a = 1, 2 
are automatically satisfied, hence by (one of the corollaries of) the Frobenius 
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theorem one can introduce curvilinear coordinates (u, v) and functions fx and 
/ 2 such that E 1 = fi(u, v)du and E 2 = f 2 (u, v)dv. Furthermore, it is possible 
to show that the functions fi and / 2 are such that fi—f 2 — A{u) + B(v), 
where A{u) and B(v) are functions of one variable, from which it follows 
that the metric can be re-written as ds 2 = (A(u) + B{v))(du 2 + dv 2 ). Note 
that the explicit forms for A(u) and B(v) follow from the vanishing of the 
curvature two-form. Then the following (differential) invariants can be used 
to solve the classification problem (for a complete solution refer to [U [38]): 
Ai(u,v) = — Tu 1 , A 2 (u,v) = — r^ 1 . In addition, we introduce in this pa- 
per the following differential invariant: As(u,v) := ^y 1 , where F = A2/A1. 
The information provided by A 3 (u,v) simplifies the classification consider- 
ably. Ultimately, we arrive at the following classification for the orthogonal 
coordinate webs (orbits) generated by the non-trivial elements of the vector 
space /C 2 (E 2 ). 



However in applications, particularly, those pertinent to the study of Hamil- 
tonian systems of classical mechanics [21 QUI EZl ESI [221 EE] , the Killing tensors 
in question arise as functions on the cotangent bundle of the configuration 
space that define principal parts of first integrals of motion. Commonly the 
components of such Killing tensors are given in terms of (cartesian) position 
coordinates. In this view it is convenient to employ a different version of the 
moving frames method [121 [13 [33] upon noticing that the Lie group SE{2) 
acts transitively on the bundle of frames. We conclude therefore that the 
same procedure can be carried out "in the group". Indeed, let x = (xi,^) 
be cartesian coordinates of E 2 . Solving the generalized Killing tensor equa- 
tion (TTj) for p = 2, n = with respect to these coordinates yields: 



where d\ = d 2 = gf- and /3i,...,/3g are arbitrary constants (of in- 
tegration) that represent the dimension of the vector space /Cq(E 2 ). The 
three-dimensional Lie group SE(2) of (orientation-preserving) isometries of 



A x = 0, A 2 = cartesian 

AxA 2 = polar 

Ai A 2 ^ 0, A 3 = parabolic 

AiA 2 ^ 0, A 3 ^ elliptic-hyperbolic 



(6) 



K 



(/3 a + 2/5 4 x 2 + fcxDdx d x 
+ ((3 2 + 2p 5 x 1 +p 6 x 2 1 )d 2 Qd 2 , 



(7) 



7 



E 2 can be represented in matrix form (using the notations adapted in |18j): 
SE(2) = |m G GL(3, R) I M = ^ * ^ ,teR 2 Re SO(2)j . In view of 
the formula <£B), the action SE(2) O E 2 given by 

x = xcosp 3 -ysinp 3 +p 1 , , , 

y = x sinp 3 + y cosp 3 + P2 

induces the corresponding action SE(2) O /Cq(E 2 ) represented by the follow- 
ing formulas [l3j [26] : 

Pi = Pi cos 2 p 3 - 2/3 3 cos p 3 sin p 3 + p 2 sin 2 p 3 - 2p 2 p A cos p 3 

-2p 2 A>sinp 3 + p 6 pl, 
p 2 = Pi sin 2 p 3 - 2p 3 cos p 3 sin p 3 + /3 2 cos 2 p 3 - 2pi/3 5 cos p 3 

+2piP 4 smp 3 + PqpI, 
P3 = (Pi - P2) sinp 3 cosp 3 + /3 3 (cos 2 p 3 - sin 2 p 3 ) (9) 

+(pi/3 4 + P2/55) cosp 3 + (pip 5 - P2P4) smp 3 - p e pip 2 , 
Pa = p A cos p 3 + /3 5 sin p 3 - p 6 p 2 , 
P_5 = P5 cos p 3 - Pa sin p 3 - p 6 pi , 
Pe = Pe- 

Using the method of moving frames "a la Fels and Olver" [T2l [131 [33] , one 
arrives at the following (algebraic) invariants (see [38] and the relevant ref- 
erences therein for more details): 

a; = p 6 , 

A' 2 = P 6 (Pi + P 2 )-Pl-Pl (10) 
A 3 = (Pe(Pi-P2)-Pi + P!) 2 + mP3 + P^) 2 - 

The fact that Ai is an invariant is obvious, the fundamental invariant A 3 
was derived for the first time in [?3] (see also [26] where A 3 was rederived 
by making use of Lie's method of infinitesimal generators). Combining the 
invariants (TTU|) . one can also make use of the invariant k 2 which is (half-) the 
distance between the singular points (foci) in the elliptic-hyperbolic case and 
is given by 

We also note that the singular points are the foci of the confocal conies (el- 
lipses and hyperbolas) that form the coordinate web in this case. In terms 
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of the algebraic invariants (TTUj) the classification §6§ can be equivalently re- 
formulated as follows [121 ES] (see also [So]). 



4 Main theorem 

The next natural step in this study is to consider the action of G on the 
product of the vector space and the underlying manifold: /Cq(M) x M or on 
n copies of the vector space: /Cq(M) x /Cq(M) x • ■ ■ x /Cq(M). In the former 
case the action leads to the covariants of Killing tensors, while in the latter - 
the joint invariants of Killing tensors (see [36J for more details). This aspect 
of the theory also has much similarity with the study of covariants and joint 
invariants in classical invariant theory (see Boutin [5J, Olver [53], Weyl [12] 
and the relevant references therein for more details). More specifically, in 
this work we extend the study of joint invariants introduced in [36]. To this 
end, we shall introduce an analogue of the concept of a resultant in classical 
invariant theory, as well as to formulate and prove an analogue of the Weyl 
theorem on joint invariants [551 |4"2"] concerning the joint action of SE(n) or 
0{n) on n copies of R n : R"xR"x-x W 1 . 

Definition 4.1 A three-dimensional orbit O G /C 2 (E 2 ) / 'SE(2) is said to be 
non- degenerate iff along O the invariant k 2 = a/A^/(A / 1 ) 2 ^ 0, where A[ 
and Ag are given by (TUP- The action SE(2) O /C 2 (E 2 ) for which k 2 ^ is 
also said to be non- degenerate. 

Recall, that all other orbits of /C 2 (E 2 ) / 'SE{2) are various degeneracies 
of the non-degenerate orbits defined above [261 EH]- The (eigenvectors of) 
Killing two-tensors corresponding to the non-degenerate orbits generate elliptic- 
hyperbolic coordinate webs [261 Ell H5J. In what follows we shall inves- 
tigate geometric meaning of the joint invariants of non-degenerate orbits 
of the orbit space (/C 2 (E 2 ) x /C 2 (E 2 ) x ••• x K? {¥?))/ SE{2). For simplic- 
ity we investigate first the joint invariants of the non-degenerate action 
SE(2) O /C 2 (E 2 ) x /C 2 (E 2 ). In this case the group SE(2) acts on each copy 
of /C 2 (E 2 ) with three-dimensional non-degenerate orbits. The group action 



a; = 0, A' 2 = 
a; = 0, A^ 

a; ^ o, a^ = o 
a;^o,a^o 



cartesian 
parabolic 
polar 

elliptic- hyperbolic 



(12) 
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is given by 




&i = 


cti cos p 3 — 2a 3 cosp 3 sinp 3 + a 2 sin p 3 — 2p 2 a 4 cosp 3 




-2p 2 a 5 smp 3 + aePa, 


d 2 = 


• 2 o ■ 2 o 

«i sm p 3 — 2« 3 cosp 3 smp 3 + a 2 cos p 3 — 2pia 5 cosp 3 




+2pia 4 sinp 3 + a 6 pi, 


d 3 = 


1 \ ■ i f 2 • 2 \ 

(ai - a 2 ) smp 3 cosp 3 + a 3 (cos z p 3 - sm p 3 ) 




+ (pi«4 + P2«s) cosp 3 + (pia s - p 2 «4) sinp 3 - a 6 pip 2 , 


ft 4 = 


0*4 COSp 3 + CK5 Sinp 3 — (^6^2; 


QJ 5 = 


Ct,5 tub — C14 bill p% — (jlqP\ j 




CK6, 


ft, — 
Pi — 


Pi cos p 3 — zp 3 cos p 3 sm p 3 -+- p 2 sm p 3 — zp 2 pA cos p 3 




-2p 2 A;smp 3 + /3 6 pi, 


A = 


fix sin 2 p 3 - 2/? 3 cosp 3 sinp 3 + f3 2 cos 2 p 3 - 2pi/5 5 cosp 3 




+2p 1/ 3 4 sinp 3 + Pepl, 


A = 


(Pi - Aj) sinp 3 cosp 3 + /3 3 (cos 2 p 3 - sin 2 p 3 ) 




+(piAi + P2A5) cosp 3 + (pips - P2/34) sinp 3 - A>pip 2 , 


h = 


Pi cosp 3 + /3 B sinp 3 - /3 6 p 2 , 


05 = 


/3 5 cosp 3 - p 4 sinp 3 - Pepi, 


h = 





(13) 



where «j, /3j, i = 1, . . . , 6 are the parameters of the respective vector spaces 
and the conditions 

2 _ a/ (oj - a| + « 6 (a 2 - «i)) 2 + 4(« 6 « 3 + a 4 a 5 ) 1 

fc2 = >/(/% - gg + - gQg + 4(ftft + /? 4/ 9 5 ) 2 ^ Q 

hold true. Six joint invariants of the action SE(2) O /C 2 (E 2 ) x /C 2 (E 2 ) follow 
immediately from (fTUj) : 

Ai = a 6 , 



5' 

2 

) 



A' 2 = a^ai + a 2) — Oi\ — a\ 

A 3 = (a 6 (ai - a 2 ) - a| + a 2 ) 2 + 4(a 6 a 3 + a 4 a 5 

A 4 = As, 

A' 5 = 

A^ = ( / 3 6 (/3i-/3 2 )-/3 4 2 + /3 2 ) 2 + 4(/3 6 /3 3 + A/? 5 ) 2 



(14) 
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The joint invariants are functionally independent (they are obtained via the 
moving frames method), any analytic function of the invariants (1141) is a 
joint invariant of the non-degenerate group action SE{2) O /Cq(E 2 ) x /Cq(E 2 ) 
(see Olver [33] for more details). In view of the Fundamental Theorem on 
invariants of regular Lie group action [33], we need to produce in total 12 (the 
dimension of the product space /Cq(E 2 ) x /Cq(E 2 )) - 3 (the dimension of the 
orbits) = 9 fundamental joint invariants. To derive the remaining three joint 
invariants, we employ geometric reasoning. Recall that each element of a 
non-degenerate orbit of the action SE{2) O /Cq(E 2 ) corresponds to a Killing 
tensor whose eigenvectors generate an elliptic-hyperbolic coordinate web. In 
turn, the coordinate web is characterized by the foci of the confocal conies 
(the two families of ellipses and hyperbolas). If F u F 2 G SE{2)/SO{2) = E 2 
are the two foci of such a coordinate web, then their respective coordinates 
(xi, i/i)f 1 and (22, D2)f 2 are given in terms of the parameters % = 1, . . . , 6 
that determine the Killing tensor via the formula ([7]) by [26J: 



Fi 



,1/2 / . \ 1/2 

-to , 1 V A 6-^ \ -Ik, 1 1 



(£2,2/2) 



F 2 



to + to 1 2 j ' to + to \ 2 



, , > 1/2 / , \ 1/2 

-to 1 / VK~^ \ 1 UfK+a^ 



(15) 



to to \ 2 ■> to to 



where o\ = /3| — /3 2 + /3 6 (/? 2 — Pi) and A 6 is given by (jl~4")) . In the case 
of the non-degenerate action SE(2) O /Cq(E 2 ) x /Cq(E 2 ), we have two more 
(generic) foci F 3 and F 4 corresponding to the second copy of /Cq(E 2 ). Their 
coordinates are given accordingly by 



F 3 



1 \ I/ 2 / / \ 1/2 

'>,, ' «6 \ 2 I ' a 6 «6 



(16) 



(£4,2/4) 



F 4 



/ \ I/ 2 / , \ 1/2 

Q6 «6 1 2 / ' «6 Q?6 
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where a 2 = ot\ — a>\ + ot & {a 2 — «i) and A' 3 is given by CHI) . Now we can treat 
the non-degenerate action SE{2) O /Cq(E 2 ) x /Cq(E 2 ), as the free and regular 
action SE{2) OE 2 xE 2 xE 2 x E 2 , where the foci F u F 2 , F 3 and F 4 belong 
to the corresponding copies of E 2 . Recall [33l 02], that any joint invariant 
J~(Fx, F 2 , -F3, F4) of the non-transitive action so defined can be written as a 
function of the interpoint distances d(Fi,F 2 ), d(F 2 ,F 3 ), d(F 3 ,F i ), di^F^Fi) 
(this result is also known as "the Weyl theorem on joint invariants" [331 
Therefore we choose the following square distances as the remaining three 
fundamental joint invariants: 

A' 7 = d 2 (F 2 , F 3 ) = (x 2 - x 3 f + (y 2 - y 3 )\ 

A> 8 = d\F,, F 3 ) = (xi - x 3 f + (yi - y 3 ) 2 , (17) 

A' 9 = d 2 (F 2 , F 4 ) = (x 2 - x 4 ) 2 + (y 2 - y 4 ) 2 , 

where (xi,yi),i = 1,2,3,4 are specified by ffTB^) and ffToT) . Direct verification 
using MAPLE shows that dA[ A dAf, A ... A dAg ^ 0, that is the nine joint 
invariants given by (fl4l) and (fT7|) are functionally independent at a generic 
point (this conclusion also follows from a geometric argument). Therefore 
we arrive at the following result. 

Theorem 4.1 Every joint invariant of the non- degenerate action SE(2) O 
/Cq(E 2 ) x/Cq(E 2 ) is a function of the nine fundamental joint invariants A-(oj , j3j) , 
i = 1, . . . , 9, j = 1, . . . , 6 given by [Lfy and ( f7?| ). 

Theorem 14.11 can naturally be extended to the general case of the non- 
degenerate group action of SE{2) on n > 2 copies of the vector space /Cq(E 2 ) 
as well as, more generally, to the case of Killing two-tensors defined on 
pseudo-Riemannian spaces of constant curvature of higher dimensions. It 
must also be mentioned that other joint invariants having a geometric mean- 
ing are the various angles, as well within the quadrilateral FiF 2 F 3 F 4: . 
For example, the angle <fi given by 

F\F 3 ■ FiF 2 

COST 



d(F 1 ,F 3 )d(F 1 ,F 2 ) 

is a joint invariant. 

We recall now the notion of a resultant in classical invariant theory (see 
Olver [33] for more details). Thus, let 

P(x) = a m x m +Ji m -ix m ~ l y H + d y m , 

Q(x) = b n x n + b n ^x n - l y + • • • + b y n [ ' 
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be two homogeneous polynomials of degrees m and n respectively. Then 
a joint invariant of the polynomials given by (j!8p is a function J(a, b) of 
the coefficients ciq, . . . , a m , b , . . . , b n preserved under the action of GL(2, M. 2 ) 
(or its subgroups). Particularly important joint invariants in this study are 
those whose vanishing is equivalent to the fact that P and Q have common 
roots. Such joint invariants are said to be resultants of the system ffTS]) . 
This important concept can be transferred naturally to the study of Killing 
tensors. 

Definition 4.2 Consider the non-degenerate action SE{2) O /Cq(E 2 ) x /Cq(E 2 ). 
Let Ki, K 2 G /Cq(E 2 ) be two Killing tensors belonging to non- degenerate or- 
bits. Then a resultant 7^[K 1 ,K 2 ] is a joint invariant of the action with the 
property that the vanishing 7?.[Ki,K 2 ] = is equivalent to the fact that the 
orthogonal coordinate webs generated by Ki and K 2 have a common focus. 

Thus, for example, the joint invariants given by (|17|) are resultants, while 
A3 + A 4 , where A 3 and A 4 are given by (fTij) . is not. Another resultant is 

A' 10 = d 2 (F u F 4 ) = { Xl - x 4 ) 2 + {yi - y A ) 2 . (19) 

In the next section we shall demonstrate how the joint invariants, in par- 
ticular, the resultants manifest themselves in the study of superintegrable 
systems. 

5 An application to the theory of superinte- 
grable systems 

In order to link the results presented in the preceding sections with the theory 
of superintegrable systems let us consider now a Hamiltonian system defined 
in E 2 by a natural Hamiltonian 

H(x,p)= 1 -( P j+ P i) + V(x) 1 (20) 

where x = (xi,x 2 ) (position coordinates) p = (pi,p 2 ) (momenta coordi- 
nates). We assume that the Hamiltonian system (T2Q|) admits two distinct in 
a certain sense (see below) first integrals Fi,F 2 , which are quadratic in the 
momenta 

F t fcp) = Kl 3 (x) PiPj + U i (x), e,i,j = 1,2 (21) 



13 



and functionally-independent with H. It is well-known that the vanishing of 
Poisson brackets {Fi, H} = {F 2 , H} = yields two sets of conditions, namely 
the Killing tensor equations 

[K 1 ,g] = [K 2) g] = 0, 

and the compatibility conditions (also known as the Bertrand-Darboux PDEs) 

d(KidV) = d(K 2 dV) = 0, 

where the components of the Killing tensors Ki, K 2 are determined by the 
quadratic in the momenta terms of the first integrals Fi and F 2 respec- 
tively given by fT2"Tj) . while the components of the (1, l)-tensors Ki, K 2 
are as follows: Kgj = Kn m g m \ I = 1,2, where are the components of 
the metric tensor that determines the kinetic part of the Hamiltonian (120]) . 
These assumptions afford orthogonal separation of variables in the associated 
Hamilton- Jacobi equation of the Hamiltonian system defined by (1201) . More 
specifically, the orthogonal coordinates are determined by the eigenvectors 
(eigenvalues) of the Killing tensors Kx and K 2 . Thus, the orthogonal separa- 
ble coordinate systems can be used to find exact solutions to the Hamiltonian 
system defined by (1201) via solving the Hamilton- Jacobi equation by separa- 
tion of variables. Moreover, in view of the above, the Hamiltonian system is 
superintegrable and multiseparable. 

Furthermore, we assume that the Killing tensors K x and K 2 belong to two 
distinct non-degenerate orbits, namely their respective eigenvectors (eigen- 
values) generate elliptic-hyperbolic coordinate webs. Let Fi, F 2 be the foci 
of the first elliptic-hyperbolic coordinate system generated by Ki, while F$, 
_F 4 - the foci of the second one, generated by K 2 . Thus, all four foci are 
distinct. We can now treat the pair {Ki,K 2 } as an element of the prod- 
uct space /Cq(E 2 ) x /Cq(E 2 ). Moreover, it is an element of a non-degenerate 
orbit of the action SE(2) O /^(E 2 ) x /C^E 2 ). Indeed, since Fx ^ F 2 and 
F3 7^ F4, the corresponding invariants k\ and k 2 given by (ITT)) do not vanish. 
Moreover, k\. In addition, we also have three joint invariants given by 
(I17p . In the most general case all of the five joint invariants kf, k 2 , A' 7 , Ag, Ag 
that completely characterize the pair {K l5 K 2 } are, in view of Theorem 14. 1[ 
functionally independent. At the same time we can treat the Killing tensors 
Ki, K 2 , as well as any linear combination C1K1 + c 2 K 2 of them as elements 
of the vector space /Cg(E 2 ). Consider now the following Killing tensor: 

K g = K x + K 2 + £g, i e K, (22) 



14 



where g is the Euclidean metric of the Hamiltonian (1201) . Clearly 

d(K ff dV) = 0, (23) 

where V is the potential part of (120]) . Our next observation is that the Killing 
tensor K 9 depends upon exactly six parameters, that is the five joint invari- 
ants and £. This is only possible if K 9 is the most general Killing tensor (j7|). 
But then the compatibility condition (1231 yields that V is constant. Since 
k\k\ 7^ 0, we conclude therefore that one of the joint invariants A' 7 , A g , A' g 
must be a vanishing resultant in order for V to be non-constant. At the same 
time two of the joint invariants cannot be resultants simultaneously, since 
k\ 7^ k\. Without loss of generality let us assume that A 8 = 7£[Ki, K 2 ] = 0. 
It follows that A' 7 = k\ and the Killing tensor K.i + K 2 in (1221) depends upon 
three parameters, namely k\ = A' 7 , k\ and Ag. We also note that since A 8 
is a vanishing resultant, the elliptic-hyperbolic coordinate webs generated 
by Ki and K 2 have a common focus. Now let us set in (j7j) the parameters 
/?3 = and f3i — fa- The resulting Killing tensor is of the form 

K' g = K[ + (3i g , (24) 

where as before g is the Euclidean metric of the Hamiltonian (120]) . Thus, the 
Killing tensor that appears in (1241) is given by 

Ki = (2{3 4 x 2 + fox\)di d x 

+(-/3 4 xi - [3 5 x 2 - /3 e xix 2 )d 1 d 2 (25) 
+(2/5 5 xx + (3 6 xl)d 2 © d 2 . 

We assume that the parameters (3^ Pb, Pe are such that both fi\ + /3| ^ 
and Pq 0. Note that /3| + P\ and Pq are invariants of the action SE{2) O 
/Cq(E 2 ) in this case. We also note that in this case the Killing tensor K' x 
generates an elliptic- hyperbolic web characterized by two foci: F\ with the 
coordinates (0,0) and F 2 with the coordinates "I") • Consider now 

another Killing tensor of the same type, namely given by 

K' 2 = (2/%z: 2 + P' 6 x 2 2 )di Qdi 

+(-^i - ^x 2 - P' & x x x 2 )di 3 2 (26) 
+(2p' 5 x 1 + P' 6 x 2 1 )d 2 ed 2 , 

under the additional assumptions that the three parameters P' A) P' 5 , P' 6 are 
such that (P' 4 ) 2 + (P' 5 ) 2 7^ and P' 6 ^ 0. Moreover, we assume in addition 
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that (/3l + f3l)//3 6 ^ [(/%) 2 + (fi' B ) 2 ]//3' e , so that the Killing tensors Ki, K' 2 ' 
are truly distinct. Note that the coordinates of the foci F 3 , F 4 of the elliptic- 

hyperbolic coordinate web generated by K' 2 are (0,0) and ( J re- 

spectively Thus, it follows that the elliptic-hyperbolic coordinate webs gen- 
erated by K' : and K 2 have one common focus at (0, 0), or, in other words, the 
pair {K'^Kj} € /Cg(E 2 ) x /Cq(E 2 ) admits a vanishing resultant. Therefore 
{K^K^} have exactly the same geometric properties as the pair of Killing 
tensors {Ki,K 2 } given by ( 12TT) . namely the elliptic-hyperbolic coordinate 
webs share one focus and the system is characterized by three parameters. 
Hence, we can identify K4 + K 2 with K' x + K' 2 . Moreover, we see that the 
sum K'j + K' 2 is given by the same formula as either of the Killing tensors 
(1251) or (1261) . In this view in order to determine the potential V in (1201) . we 
substitute the formula (1251) into the compatibility condition (1231) and solve 
the resulting PDE for V. More specifically, after the substitution we set suc- 
cessively in ( 123]) /3 4 = /? 5 = 0, which defines a polar web with singular point 
at (0, 0), then /3 4 = /5 6 = 0, which defines a parabolic web with singular point 
at (0, 0) and the focal axis, and finally (3$ = (3q = 0, which defines 

a parabolic web with singular point at (0, 0) and focal axis the ?/-axis. This 
results in three PDEs whose unique solution can be easily found. Indeed, the 
PDEs enjoy the following forms: 

2x 2 V Xl - 2x{V X2 - xix 2 {V X2X2 - V XlXl ) + (x\ - x\)V XlX2 = 0, (27) 

3V X2 + x 2 (V X2X2 - V X1X1 ) + 2 Xl V XlX2 = 0, (28) 

3V Xl + Xl {V X2X2 - V X1X1 ) + 2x 2 V XlX2 = (29) 
The PDEs (|2Tl)-(p9l) yield 

x 2 V Xl - XiV X2 = 0, (30) 

which implies that (in terms of polar coordinates) 

V(r, 9) = V(r cos 9, r sin 9) (31) 

is independent of 9. Now the general solution of (!29l in polar coordinates 
has the form 

V{r,6) = F(r) + (32) 
where F and G are arbitrary functions. Thus by (130!) we have 

G{9) = i (33) 
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where t is an arbitrary constant. Substituting (1521) into (13"U1) (in polar coor- 
dinates) and taking (1331) into account, we find that 



F(r) = ---- + n, (34) 

where £ as in (133]) and m, n are arbitrary constants. Substituting (133]) and 
(134]) into (1321) and transforming back to cartesian coordinates, we arrive at 
the potential of the Kepler problem: 

V(xi,x 2 ) = , , (35) 

where without loss of generality we set m = —1, n = 0. The inverse problem, 
namely when one starts with the Kepler potential ( 1351) and then finds the 
most general Killing tensor compatible with (135!) via (|23|) . can be solved in 
the same manner, that is by solving the corresponding (Bertrand-Darboux) 
PDE. The calculations are straightforward and we present here the result 
only. Thus, the most general Killing tensor compatible with the Kepler 
potential ( 135]) is precisely the three-parameter family of Killing tensors given 
by (l25|) . We conclude therefore that we have proven the following theorem. 



Theorem 5.1 Let the potential V of the general Hamiltonian be com- 
patible via (d3|) with any two non- degenerate Killing tensors K 1; K 2 G /Cq(E 2 ). 
Then the following statements are equivalent. 

(1) The pair of Killing tensors {K 1; K 2 } G /Cq(E 2 ) X /Cq(IE 2 ) admits one 
vanishing resultant 72. [Ki, K 2 ] . 

(2) The potential V given by l[2lfy is the Kepler potential (E^). 



6 Conclusions 

The results presented in this article can naturally be generalized by increasing 
the dimension of the underlying (pseudo-) Riemannian manifold, changing 
its curvature and the signature of its metric. Some of these cases will be 
investigated by the authors in the forthcoming papers on the subject. 
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